*
* $Id$
*
* $Log: gphrin.F,v $
* Revision 1.1.1.1  2002/07/24 15:56:26  rdm
* initial import into CVS
*
* Revision 1.1.1.1  2002/06/16 15:18:41  hristov
* Separate distribution  of Geant3
*
* Revision 1.1.1.1  1999/05/18 15:55:20  fca
* AliRoot sources
*
* Revision 1.1.1.1  1995/10/24 10:21:29  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.22  by  S.Giani
*-- Author :
      SUBROUTINE G3PHRIN
C.
C.    ******************************************************************
C.    *                                                                *
C     *       Initialize decay mode rates' constants for               *
C.    *       the photoelectric effect                                 *
C.    *                                                                *
C.    *                                                                *
C.    *    ==>Called by : G3PHINI                                      *
C.    *       Author    J. Chwastowski                                 *
C.    *                                                                *
C.    ******************************************************************
C.
      IMPLICIT NONE
#include "geant321/gcpmxz.inc"
#include "geant321/gcphpr.inc"
#include "geant321/gcphnr.inc"
#include "geant321/gcphrd.inc"
      INTEGER I,J,ITRN
      REAL TMP,TMP1
      DIMENSION TMP(ISHLS,MAXELZ),TMP1(4,MAXELZ)
      DIMENSION ITRN(29)
      DATA ITRN / 3,4,6, 7,11,12,
     +            3,4,6, 7,11,12,18,19,
     +            4,5,8,10,13,17,20,
     +            5,8,9,10,13,14,17,20 /
      DATA (TMP(J,1),J=1,ISHLS) / ISHLS*0.0 /
      DATA (TMP(J,2),J=1,ISHLS) / ISHLS*0.0 /
      DATA (TMP(J,3),J=1,ISHLS) / ISHLS*0.0 /
      DATA (TMP(J,4),J=1,ISHLS) / ISHLS*0.0 /
      DATA (TMP(J,  5),J=1,ISHLS) /
     +  0.333E+00, 0.667E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.333E+00, 0.667E+00
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
      DATA (TMP(J,  6),J=1,ISHLS) /
     +  0.334E+00, 0.666E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.332E+00, 0.668E+00
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
      DATA (TMP(J,  7),J=1,ISHLS) /
     +  0.334E+00, 0.666E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.332E+00, 0.668E+00
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
      DATA (TMP(J,  8),J=1,ISHLS) /
     +  0.334E+00, 0.666E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.331E+00, 0.669E+00
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
      DATA (TMP(J,  9),J=1,ISHLS) /
     +  0.334E+00, 0.666E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.331E+00, 0.669E+00
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
      DATA (TMP(J, 10),J=1,ISHLS) /
     +  0.335E+00, 0.665E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.329E+00, 0.671E+00
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
      DATA (TMP(J, 11),J=1,ISHLS) /
     +  0.335E+00, 0.665E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.328E+00, 0.672E+00
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.233E-09, 0.100E+01
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.100E+01,  7*0.0 /
      DATA (TMP(J, 12),J=1,ISHLS) /
     +  0.335E+00, 0.665E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.326E+00, 0.674E+00
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.159E-09, 0.100E+01
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.100E+01,  7*0.0 /
      DATA (TMP(J, 13),J=1,ISHLS) /
     +  0.333E+00, 0.661E+00, 0.203E-02, 0.404E-02
     +, 0.000E+00, 0.000E+00, 0.291E+00, 0.608E+00
     +, 0.342E-01, 0.663E-01, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.239E-09, 0.100E+01
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.100E+01,  7*0.0 /
      DATA (TMP(J, 14),J=1,ISHLS) /
     +  0.330E+00, 0.654E+00, 0.562E-02, 0.111E-01
     +, 0.000E+00, 0.000E+00, 0.216E+00, 0.456E+00
     +, 0.111E+00, 0.216E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.408E-09, 0.100E+01
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.100E+01,  7*0.0 /
      DATA (TMP(J, 15),J=1,ISHLS) /
     +  0.325E+00, 0.645E+00, 0.100E-01, 0.199E-01
     +, 0.000E+00, 0.000E+00, 0.141E+00, 0.300E+00
     +, 0.190E+00, 0.369E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.703E-09, 0.100E+01
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.100E+01,  7*0.0 /
      DATA (TMP(J, 16),J=1,ISHLS) /
     +  0.321E+00, 0.635E+00, 0.149E-01, 0.296E-01
     +, 0.000E+00, 0.000E+00, 0.881E-01, 0.191E+00
     +, 0.246E+00, 0.475E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.120E-08, 0.100E+01
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.100E+01,  7*0.0 /
      DATA (TMP(J, 17),J=1,ISHLS) /
     +  0.315E+00, 0.625E+00, 0.201E-01, 0.398E-01
     +, 0.000E+00, 0.000E+00, 0.559E-01, 0.124E+00
     +, 0.281E+00, 0.540E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.199E-08, 0.100E+01
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.100E+01,  7*0.0 /
      DATA (TMP(J, 18),J=1,ISHLS) /
     +  0.314E+00, 0.611E+00, 0.253E-01, 0.500E-01
     +, 0.000E+00, 0.000E+00, 0.368E-01, 0.828E-01
     +, 0.302E+00, 0.579E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.324E-08, 0.100E+01
     +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.100E+01,  7*0.0 /
      DATA (TMP(J, 19),J=1,ISHLS) /
     +  0.307E+00, 0.606E+00, 0.293E-01, 0.580E-01
     +, 0.000E+00, 0.000E+00, 0.260E-01, 0.597E-01
     +, 0.313E+00, 0.601E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.480E-08, 0.976E+00
     +, 0.000E+00, 0.238E-01, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.976E+00, 0.000E+00, 0.000E+00
     +, 0.238E-01,  4*0.0 /
      DATA (TMP(J, 20),J=1,ISHLS) /
     +  0.304E+00, 0.600E+00, 0.323E-01, 0.639E-01
     +, 0.000E+00, 0.000E+00, 0.193E-01, 0.455E-01
     +, 0.321E+00, 0.614E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.678E-08, 0.929E+00
     +, 0.000E+00, 0.709E-01, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.929E+00, 0.000E+00, 0.000E+00
     +, 0.709E-01,  4*0.0 /
      DATA (TMP(J, 21),J=1,ISHLS) /
     +  0.303E+00, 0.597E+00, 0.335E-01, 0.661E-01
     +, 0.000E+00, 0.000E+00, 0.155E-01, 0.376E-01
     +, 0.326E+00, 0.621E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.504E-08, 0.452E+00
     +, 0.513E+00, 0.347E-01, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.461E+00, 0.506E-01, 0.453E+00
     +, 0.353E-01,  4*0.0 /
      DATA (TMP(J, 22),J=1,ISHLS) /
     +  0.303E+00, 0.595E+00, 0.343E-01, 0.677E-01
     +, 0.000E+00, 0.000E+00, 0.128E-01, 0.317E-01
     +, 0.330E+00, 0.626E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.463E-08, 0.276E+00
     +, 0.703E+00, 0.208E-01, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.284E+00, 0.699E-01, 0.624E+00
     +, 0.214E-01,  4*0.0 /
      DATA (TMP(J, 23),J=1,ISHLS) /
     +  0.302E+00, 0.594E+00, 0.350E-01, 0.690E-01
     +, 0.000E+00, 0.000E+00, 0.107E-01, 0.273E-01
     +, 0.333E+00, 0.629E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.471E-08, 0.191E+00
     +, 0.795E+00, 0.139E-01, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.198E+00, 0.793E-01, 0.708E+00
     +, 0.145E-01,  4*0.0 /
      DATA (TMP(J, 24),J=1,ISHLS) /
     +  0.303E+00, 0.594E+00, 0.348E-01, 0.685E-01
     +, 0.000E+00, 0.000E+00, 0.926E-02, 0.245E-01
     +, 0.335E+00, 0.631E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.475E-08, 0.130E+00
     +, 0.861E+00, 0.888E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.136E+00, 0.863E-01, 0.769E+00
     +, 0.910E-02,  4*0.0 /
      DATA (TMP(J, 25),J=1,ISHLS) /
     +  0.302E+00, 0.592E+00, 0.359E-01, 0.708E-01
     +, 0.000E+00, 0.000E+00, 0.778E-02, 0.213E-01
     +, 0.337E+00, 0.633E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.573E-08, 0.111E+00
     +, 0.881E+00, 0.759E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.117E+00, 0.882E-01, 0.787E+00
     +, 0.803E-02,  4*0.0 /
      DATA (TMP(J, 26),J=1,ISHLS) /
     +  0.302E+00, 0.591E+00, 0.363E-01, 0.714E-01
     +, 0.000E+00, 0.000E+00, 0.675E-02, 0.192E-01
     +, 0.339E+00, 0.635E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.657E-08, 0.900E-01
     +, 0.904E+00, 0.597E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.955E-01, 0.906E-01, 0.807E+00
     +, 0.634E-02,  4*0.0 /
      DATA (TMP(J, 27),J=1,ISHLS) /
     +  0.302E+00, 0.590E+00, 0.366E-01, 0.719E-01
     +, 0.000E+00, 0.000E+00, 0.589E-02, 0.176E-01
     +, 0.341E+00, 0.635E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.767E-08, 0.751E-01
     +, 0.920E+00, 0.483E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.803E-01, 0.924E-01, 0.822E+00
     +, 0.515E-02,  4*0.0 /
      DATA (TMP(J, 28),J=1,ISHLS) /
     +  0.301E+00, 0.589E+00, 0.369E-01, 0.724E-01
     +, 0.000E+00, 0.000E+00, 0.521E-02, 0.162E-01
     +, 0.343E+00, 0.635E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.902E-08, 0.640E-01
     +, 0.932E+00, 0.398E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.690E-01, 0.938E-01, 0.833E+00
     +, 0.428E-02,  4*0.0 /
      DATA (TMP(J, 29),J=1,ISHLS) /
     +  0.302E+00, 0.590E+00, 0.366E-01, 0.717E-01
     +, 0.000E+00, 0.000E+00, 0.470E-02, 0.153E-01
     +, 0.345E+00, 0.635E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.104E-07, 0.537E-01
     +, 0.943E+00, 0.325E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.585E-01, 0.954E-01, 0.845E+00
     +, 0.134E-02,  4*0.0 /
      DATA (TMP(J, 30),J=1,ISHLS) /
     +  0.302E+00, 0.588E+00, 0.373E-01, 0.730E-01
     +, 0.000E+00, 0.000E+00, 0.416E-02, 0.141E-01
     +, 0.347E+00, 0.635E+00, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.000E+00, 0.127E-07, 0.488E-01
     +, 0.948E+00, 0.285E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.532E-01, 0.956E-01, 0.848E+00
     +, 0.311E-02,  4*0.0 /
      DATA (TMP(J, 31),J=1,ISHLS) /
     +  0.301E+00, 0.586E+00, 0.381E-01, 0.744E-01
     +, 0.217E-03, 0.414E-03, 0.365E-02, 0.131E-01
     +, 0.346E+00, 0.631E+00, 0.208E-02, 0.373E-02
     +, 0.000E+00, 0.000E+00, 0.153E-07, 0.449E-01
     +, 0.952E+00, 0.345E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.494E-01, 0.958E-01, 0.851E+00
     +, 0.378E-02,  4*0.0 /
      DATA (TMP(J, 32),J=1,ISHLS) /
     +  0.301E+00, 0.583E+00, 0.388E-01, 0.759E-01
     +, 0.633E-03, 0.122E-02, 0.322E-02, 0.121E-01
     +, 0.344E+00, 0.624E+00, 0.595E-02, 0.106E-01
     +, 0.000E+00, 0.000E+00, 0.187E-07, 0.421E-01
     +, 0.954E+00, 0.376E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.465E-01, 0.961E-01, 0.853E+00
     +, 0.414E-02,  4*0.0 /
      DATA (TMP(J, 33),J=1,ISHLS) /
     +  0.299E+00, 0.580E+00, 0.396E-01, 0.772E-01
     +, 0.121E-02, 0.232E-02, 0.282E-02, 0.112E-01
     +, 0.341E+00, 0.614E+00, 0.110E-01, 0.197E-01
     +, 0.000E+00, 0.000E+00, 0.229E-07, 0.398E-01
     +, 0.956E+00, 0.398E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.442E-01, 0.962E-01, 0.855E+00
     +, 0.441E-02,  4*0.0 /
      DATA (TMP(J, 34),J=1,ISHLS) /
     +  0.298E+00, 0.577E+00, 0.403E-01, 0.787E-01
     +, 0.191E-02, 0.368E-02, 0.249E-02, 0.105E-01
     +, 0.336E+00, 0.604E+00, 0.170E-01, 0.302E-01
     +, 0.000E+00, 0.000E+00, 0.279E-07, 0.379E-01
     +, 0.958E+00, 0.414E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.425E-01, 0.964E-01, 0.857E+00
     +, 0.462E-02,  4*0.0 /
      DATA (TMP(J, 35),J=1,ISHLS) /
     +  0.297E+00, 0.574E+00, 0.410E-01, 0.799E-01
     +, 0.273E-02, 0.526E-02, 0.220E-02, 0.979E-02
     +, 0.331E+00, 0.592E+00, 0.235E-01, 0.417E-01
     +, 0.000E+00, 0.000E+00, 0.340E-07, 0.364E-01
     +, 0.959E+00, 0.427E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.411E-01, 0.965E-01, 0.858E+00
     +, 0.480E-02,  4*0.0 /
      DATA (TMP(J, 36),J=1,ISHLS) /
     +  0.296E+00, 0.571E+00, 0.416E-01, 0.812E-01
     +, 0.366E-02, 0.704E-02, 0.194E-02, 0.921E-02
     +, 0.325E+00, 0.579E+00, 0.305E-01, 0.539E-01
     +, 0.000E+00, 0.000E+00, 0.413E-07, 0.352E-01
     +, 0.960E+00, 0.438E-02, 0.000E+00, 0.000E+00
     +, 0.000E+00, 0.400E-01, 0.966E-01, 0.858E+00
     +, 0.495E-02,  4*0.0 /
      DATA (TMP(J, 37),J=1,ISHLS) /
     +  0.295E+00, 0.567E+00, 0.423E-01, 0.824E-01
     +, 0.440E-02, 0.855E-02, 0.173E-02, 0.873E-02
     +, 0.321E+00, 0.569E+00, 0.358E-01, 0.636E-01
     +, 0.000E+00, 0.000E+00, 0.502E-07, 0.341E-01
     +, 0.961E+00, 0.463E-02, 0.000E+00, 0.130E-03
     +, 0.000E+00, 0.390E-01, 0.966E-01, 0.859E+00
     +, 0.529E-02, 0.000E+00, 0.000E+00, 0.148E-03, 0.0 /
      DATA (TMP(J, 38),J=1,ISHLS) /
     +  0.294E+00, 0.565E+00, 0.429E-01, 0.837E-01
     +, 0.506E-02, 0.985E-02, 0.155E-02, 0.836E-02
     +, 0.318E+00, 0.560E+00, 0.402E-01, 0.712E-01
     +, 0.000E+00, 0.000E+00, 0.607E-07, 0.331E-01
     +, 0.962E+00, 0.489E-02, 0.000E+00, 0.421E-03
     +, 0.000E+00, 0.382E-01, 0.967E-01, 0.859E+00
     +, 0.561E-02, 0.000E+00, 0.000E+00, 0.484E-03, 0.0 /
      DATA (TMP(J, 39),J=1,ISHLS) /
     +  0.293E+00, 0.562E+00, 0.436E-01, 0.848E-01
     +, 0.554E-02, 0.108E-01, 0.140E-02, 0.807E-02
     +, 0.317E+00, 0.554E+00, 0.432E-01, 0.762E-01
     +, 0.000E+00, 0.000E+00, 0.730E-07, 0.322E-01
     +, 0.957E+00, 0.501E-02, 0.496E-02, 0.465E-03
     +, 0.000E+00, 0.374E-01, 0.962E-01, 0.855E+00
     +, 0.579E-02, 0.494E-03, 0.437E-02, 0.535E-03, 0.0 /
      DATA (TMP(J, 40),J=1,ISHLS) /
     +  0.293E+00, 0.560E+00, 0.441E-01, 0.860E-01
     +, 0.595E-02, 0.116E-01, 0.127E-02, 0.785E-02
     +, 0.316E+00, 0.549E+00, 0.457E-01, 0.802E-01
     +, 0.000E+00, 0.000E+00, 0.871E-07, 0.312E-01
     +, 0.952E+00, 0.508E-02, 0.117E-01, 0.480E-03
     +, 0.000E+00, 0.367E-01, 0.957E-01, 0.850E+00
     +, 0.593E-02, 0.116E-02, 0.103E-01, 0.559E-03, 0.0 /
      DATA (TMP(J, 41),J=1,ISHLS) /
     +  0.292E+00, 0.558E+00, 0.447E-01, 0.871E-01
     +, 0.623E-02, 0.121E-01, 0.116E-02, 0.770E-02
     +, 0.316E+00, 0.546E+00, 0.472E-01, 0.823E-01
     +, 0.000E+00, 0.000E+00, 0.103E-06, 0.303E-01
     +, 0.941E+00, 0.506E-02, 0.232E-01, 0.198E-03
     +, 0.000E+00, 0.359E-01, 0.946E-01, 0.841E+00
     +, 0.595E-02, 0.232E-02, 0.204E-01, 0.233E-03, 0.0 /
      DATA (TMP(J, 42),J=1,ISHLS) /
     +  0.292E+00, 0.556E+00, 0.453E-01, 0.881E-01
     +, 0.656E-02, 0.127E-01, 0.107E-02, 0.758E-02
     +, 0.315E+00, 0.542E+00, 0.491E-01, 0.853E-01
     +, 0.000E+00, 0.000E+00, 0.123E-06, 0.295E-01
     +, 0.933E+00, 0.510E-02, 0.321E-01, 0.194E-03
     +, 0.000E+00, 0.352E-01, 0.939E-01, 0.833E+00
     +, 0.604E-02, 0.319E-02, 0.282E-01, 0.231E-03, 0.0 /
      DATA (TMP(J, 43),J=1,ISHLS) /
     +  0.291E+00, 0.554E+00, 0.458E-01, 0.890E-01
     +, 0.688E-02, 0.133E-01, 0.978E-03, 0.749E-02
     +, 0.315E+00, 0.538E+00, 0.508E-01, 0.878E-01
     +, 0.000E+00, 0.000E+00, 0.145E-06, 0.288E-01
     +, 0.924E+00, 0.510E-02, 0.416E-01, 0.189E-03
     +, 0.000E+00, 0.347E-01, 0.930E-01, 0.825E+00
     +, 0.610E-02, 0.413E-02, 0.364E-01, 0.226E-03, 0.0 /
      DATA (TMP(J, 44),J=1,ISHLS) /
     +  0.291E+00, 0.552E+00, 0.463E-01, 0.901E-01
     +, 0.716E-02, 0.139E-01, 0.902E-03, 0.745E-02
     +, 0.315E+00, 0.535E+00, 0.524E-01, 0.899E-01
     +, 0.000E+00, 0.000E+00, 0.172E-06, 0.281E-01
     +, 0.915E+00, 0.509E-02, 0.515E-01, 0.184E-03
     +, 0.000E+00, 0.341E-01, 0.921E-01, 0.817E+00
     +, 0.614E-02, 0.511E-02, 0.450E-01, 0.222E-03, 0.0 /
      DATA (TMP(J, 45),J=1,ISHLS) /
     +  0.291E+00, 0.550E+00, 0.468E-01, 0.909E-01
     +, 0.742E-02, 0.144E-01, 0.834E-03, 0.744E-02
     +, 0.315E+00, 0.531E+00, 0.538E-01, 0.918E-01
     +, 0.000E+00, 0.000E+00, 0.202E-06, 0.274E-01
     +, 0.906E+00, 0.508E-02, 0.617E-01, 0.177E-03
     +, 0.000E+00, 0.336E-01, 0.911E-01, 0.809E+00
     +, 0.618E-02, 0.612E-02, 0.540E-01, 0.216E-03, 0.0 /
      DATA (TMP(J, 46),J=1,ISHLS) /
     +  0.290E+00, 0.548E+00, 0.472E-01, 0.918E-01
     +, 0.762E-02, 0.147E-01, 0.773E-03, 0.747E-02
     +, 0.315E+00, 0.529E+00, 0.547E-01, 0.925E-01
     +, 0.000E+00, 0.000E+00, 0.236E-06, 0.267E-01
     +, 0.894E+00, 0.502E-02, 0.744E-01, 0.000E+00
     +, 0.000E+00, 0.332E-01, 0.900E-01, 0.799E+00
     +, 0.617E-02, 0.737E-02, 0.646E-01,  2*0.0 /
      DATA (TMP(J, 47),J=1,ISHLS) /
     +  0.290E+00, 0.547E+00, 0.477E-01, 0.927E-01
     +, 0.790E-02, 0.153E-01, 0.719E-03, 0.751E-02
     +, 0.315E+00, 0.525E+00, 0.564E-01, 0.947E-01
     +, 0.000E+00, 0.000E+00, 0.276E-06, 0.262E-01
     +, 0.886E+00, 0.503E-02, 0.828E-01, 0.164E-03
     +, 0.000E+00, 0.329E-01, 0.891E-01, 0.791E+00
     +, 0.624E-02, 0.820E-02, 0.723E-01, 0.202E-03, 0.0 /
      DATA (TMP(J, 48),J=1,ISHLS) /
     +  0.290E+00, 0.545E+00, 0.481E-01, 0.934E-01
     +, 0.823E-02, 0.159E-01, 0.668E-03, 0.758E-02
     +, 0.315E+00, 0.521E+00, 0.582E-01, 0.971E-01
     +, 0.000E+00, 0.000E+00, 0.323E-06, 0.257E-01
     +, 0.878E+00, 0.505E-02, 0.905E-01, 0.400E-03
     +, 0.000E+00, 0.326E-01, 0.884E-01, 0.784E+00
     +, 0.632E-02, 0.894E-02, 0.790E-01, 0.500E-03, 0.0 /
      DATA (TMP(J, 49),J=1,ISHLS) /
     +  0.290E+00, 0.543E+00, 0.484E-01, 0.941E-01
     +, 0.853E-02, 0.165E-01, 0.623E-03, 0.766E-02
     +, 0.315E+00, 0.516E+00, 0.598E-01, 0.995E-01
     +, 0.507E-03, 0.778E-03, 0.378E-06, 0.253E-01
     +, 0.871E+00, 0.507E-02, 0.978E-01, 0.503E-03
     +, 0.000E+00, 0.324E-01, 0.877E-01, 0.778E+00
     +, 0.641E-02, 0.965E-02, 0.855E-01, 0.634E-03, 0.0 /
      DATA (TMP(J, 50),J=1,ISHLS) /
     +  0.289E+00, 0.541E+00, 0.489E-01, 0.948E-01
     +, 0.884E-02, 0.171E-01, 0.580E-03, 0.776E-02
     +, 0.314E+00, 0.511E+00, 0.614E-01, 0.101E+00
     +, 0.140E-02, 0.218E-02, 0.439E-06, 0.249E-01
     +, 0.865E+00, 0.510E-02, 0.104E+00, 0.574E-03
     +, 0.000E+00, 0.322E-01, 0.870E-01, 0.772E+00
     +, 0.651E-02, 0.103E-01, 0.914E-01, 0.732E-03, 0.0 /
      DATA (TMP(J, 51),J=1,ISHLS) /
     +  0.289E+00, 0.540E+00, 0.492E-01, 0.953E-01
     +, 0.915E-02, 0.177E-01, 0.541E-03, 0.787E-02
     +, 0.313E+00, 0.505E+00, 0.629E-01, 0.103E+00
     +, 0.260E-02, 0.406E-02, 0.511E-06, 0.245E-01
     +, 0.859E+00, 0.513E-02, 0.111E+00, 0.634E-03
     +, 0.000E+00, 0.321E-01, 0.864E-01, 0.766E+00
     +, 0.661E-02, 0.109E-01, 0.967E-01, 0.816E-03, 0.0 /
      DATA (TMP(J, 52),J=1,ISHLS) /
     +  0.289E+00, 0.537E+00, 0.496E-01, 0.961E-01
     +, 0.945E-02, 0.183E-01, 0.505E-03, 0.800E-02
     +, 0.312E+00, 0.500E+00, 0.642E-01, 0.105E+00
     +, 0.404E-02, 0.631E-02, 0.593E-06, 0.242E-01
     +, 0.854E+00, 0.516E-02, 0.116E+00, 0.688E-03
     +, 0.000E+00, 0.320E-01, 0.859E-01, 0.762E+00
     +, 0.672E-02, 0.114E-01, 0.102E+00, 0.894E-03, 0.0 /
      DATA (TMP(J, 53),J=1,ISHLS) /
     +  0.289E+00, 0.536E+00, 0.499E-01, 0.967E-01
     +, 0.974E-02, 0.189E-01, 0.473E-03, 0.816E-02
     +, 0.311E+00, 0.494E+00, 0.656E-01, 0.107E+00
     +, 0.571E-02, 0.891E-02, 0.688E-06, 0.240E-01
     +, 0.849E+00, 0.519E-02, 0.121E+00, 0.738E-03
     +, 0.000E+00, 0.319E-01, 0.854E-01, 0.757E+00
     +, 0.683E-02, 0.119E-01, 0.106E+00, 0.968E-03, 0.0 /
      DATA (TMP(J, 54),J=1,ISHLS) /
     +  0.289E+00, 0.534E+00, 0.502E-01, 0.973E-01
     +, 0.100E-01, 0.195E-01, 0.442E-03, 0.833E-02
     +, 0.309E+00, 0.488E+00, 0.667E-01, 0.108E+00
     +, 0.759E-02, 0.118E-01, 0.792E-06, 0.237E-01
     +, 0.844E+00, 0.523E-02, 0.126E+00, 0.784E-03
     +, 0.000E+00, 0.319E-01, 0.849E-01, 0.753E+00
     +, 0.694E-02, 0.123E-01, 0.110E+00, 0.104E-02, 0.0 /
      DATA (TMP(J, 55),J=1,ISHLS) /
     +  0.288E+00, 0.533E+00, 0.505E-01, 0.978E-01
     +, 0.103E-01, 0.200E-01, 0.415E-03, 0.853E-02
     +, 0.309E+00, 0.482E+00, 0.679E-01, 0.109E+00
     +, 0.900E-02, 0.142E-01, 0.915E-06, 0.235E-01
     +, 0.840E+00, 0.526E-02, 0.130E+00, 0.857E-03
     +, 0.000E+00, 0.319E-01, 0.845E-01, 0.749E+00
     +, 0.706E-02, 0.128E-01, 0.114E+00, 0.115E-02, 0.0 /
      DATA (TMP(J, 56),J=1,ISHLS) /
     +  0.289E+00, 0.531E+00, 0.508E-01, 0.984E-01
     +, 0.106E-01, 0.206E-01, 0.391E-03, 0.876E-02
     +, 0.308E+00, 0.477E+00, 0.692E-01, 0.110E+00
     +, 0.103E-01, 0.162E-01, 0.105E-05, 0.232E-01
     +, 0.836E+00, 0.529E-02, 0.135E+00, 0.928E-03
     +, 0.000E+00, 0.320E-01, 0.841E-01, 0.745E+00
     +, 0.718E-02, 0.132E-01, 0.118E+00, 0.126E-02, 0.0 /
      DATA (TMP(J, 57),J=1,ISHLS) /
     +  0.289E+00, 0.530E+00, 0.510E-01, 0.989E-01
     +, 0.108E-01, 0.211E-01, 0.368E-03, 0.903E-02
     +, 0.308E+00, 0.472E+00, 0.705E-01, 0.112E+00
     +, 0.112E-01, 0.175E-01, 0.121E-05, 0.229E-01
     +, 0.831E+00, 0.531E-02, 0.139E+00, 0.980E-03
     +, 0.108E-02, 0.320E-01, 0.836E-01, 0.740E+00
     +, 0.730E-02, 0.135E-01, 0.121E+00, 0.134E-02, 0.103E-02 /
      DATA (TMP(J, 58),J=1,ISHLS) /
     +  0.288E+00, 0.529E+00, 0.514E-01, 0.995E-01
     +, 0.109E-01, 0.212E-01, 0.348E-03, 0.938E-02
     +, 0.310E+00, 0.471E+00, 0.710E-01, 0.111E+00
     +, 0.104E-01, 0.161E-01, 0.140E-05, 0.229E-01
     +, 0.833E+00, 0.530E-02, 0.138E+00, 0.927E-03
     +, 0.000E+00, 0.323E-01, 0.839E-01, 0.742E+00
     +, 0.736E-02, 0.135E-01, 0.120E+00, 0.128E-02, 0.0 /
      DATA (TMP(J, 59),J=1,ISHLS) /
     +  0.289E+00, 0.527E+00, 0.517E-01, 0.100E+00
     +, 0.110E-01, 0.214E-01, 0.330E-03, 0.973E-02
     +, 0.311E+00, 0.469E+00, 0.718E-01, 0.112E+00
     +, 0.105E-01, 0.160E-01, 0.161E-05, 0.228E-01
     +, 0.831E+00, 0.532E-02, 0.140E+00, 0.922E-03
     +, 0.000E+00, 0.325E-01, 0.838E-01, 0.740E+00
     +, 0.745E-02, 0.136E-01, 0.121E+00, 0.129E-02, 0.0 /
      DATA (TMP(J, 60),J=1,ISHLS) /
     +  0.289E+00, 0.526E+00, 0.520E-01, 0.101E+00
     +, 0.112E-01, 0.217E-01, 0.313E-03, 0.101E-01
     +, 0.313E+00, 0.466E+00, 0.725E-01, 0.112E+00
     +, 0.105E-01, 0.158E-01, 0.185E-05, 0.227E-01
     +, 0.830E+00, 0.532E-02, 0.141E+00, 0.918E-03
     +, 0.000E+00, 0.327E-01, 0.834E-01, 0.740E+00
     +, 0.753E-02, 0.136E-01, 0.122E+00, 0.129E-02, 0.0 /
      DATA (TMP(J, 61),J=1,ISHLS) /
     +  0.289E+00, 0.524E+00, 0.522E-01, 0.101E+00
     +, 0.113E-01, 0.219E-01, 0.297E-03, 0.105E-01
     +, 0.314E+00, 0.463E+00, 0.732E-01, 0.112E+00
     +, 0.105E-01, 0.156E-01, 0.213E-05, 0.226E-01
     +, 0.830E+00, 0.532E-02, 0.142E+00, 0.912E-03
     +, 0.000E+00, 0.330E-01, 0.834E-01, 0.739E+00
     +, 0.761E-02, 0.137E-01, 0.122E+00, 0.130E-02, 0.0 /
      DATA (TMP(J, 62),J=1,ISHLS) /
     +  0.289E+00, 0.523E+00, 0.525E-01, 0.101E+00
     +, 0.114E-01, 0.221E-01, 0.283E-03, 0.110E-01
     +, 0.316E+00, 0.461E+00, 0.740E-01, 0.112E+00
     +, 0.105E-01, 0.155E-01, 0.243E-05, 0.225E-01
     +, 0.829E+00, 0.531E-02, 0.142E+00, 0.905E-03
     +, 0.000E+00, 0.333E-01, 0.835E-01, 0.737E+00
     +, 0.772E-02, 0.138E-01, 0.123E+00, 0.131E-02, 0.0 /
      DATA (TMP(J, 63),J=1,ISHLS) /
     +  0.289E+00, 0.522E+00, 0.527E-01, 0.102E+00
     +, 0.115E-01, 0.223E-01, 0.269E-03, 0.115E-01
     +, 0.318E+00, 0.458E+00, 0.747E-01, 0.112E+00
     +, 0.104E-01, 0.153E-01, 0.277E-05, 0.224E-01
     +, 0.828E+00, 0.532E-02, 0.143E+00, 0.898E-03
     +, 0.000E+00, 0.335E-01, 0.832E-01, 0.737E+00
     +, 0.779E-02, 0.138E-01, 0.123E+00, 0.131E-02, 0.0 /
      DATA (TMP(J, 64),J=1,ISHLS) /
     +  0.290E+00, 0.521E+00, 0.529E-01, 0.102E+00
     +, 0.116E-01, 0.226E-01, 0.257E-03, 0.120E-01
     +, 0.318E+00, 0.454E+00, 0.755E-01, 0.112E+00
     +, 0.111E-01, 0.161E-01, 0.315E-05, 0.223E-01
     +, 0.826E+00, 0.533E-02, 0.145E+00, 0.933E-03
     +, 0.865E-03, 0.338E-01, 0.830E-01, 0.734E+00
     +, 0.790E-02, 0.139E-01, 0.125E+00, 0.138E-02, 0.798E-03 /
      DATA (TMP(J, 65),J=1,ISHLS) /
     +  0.290E+00, 0.520E+00, 0.532E-01, 0.103E+00
     +, 0.117E-01, 0.226E-01, 0.245E-03, 0.126E-01
     +, 0.321E+00, 0.453E+00, 0.760E-01, 0.112E+00
     +, 0.104E-01, 0.148E-01, 0.360E-05, 0.223E-01
     +, 0.827E+00, 0.532E-02, 0.144E+00, 0.887E-03
     +, 0.000E+00, 0.342E-01, 0.832E-01, 0.735E+00
     +, 0.799E-02, 0.139E-01, 0.124E+00, 0.133E-02, 0.0 /
      DATA (TMP(J, 66),J=1,ISHLS) /
     +  0.290E+00, 0.519E+00, 0.534E-01, 0.103E+00
     +, 0.118E-01, 0.228E-01, 0.234E-03, 0.133E-01
     +, 0.323E+00, 0.450E+00, 0.768E-01, 0.112E+00
     +, 0.104E-01, 0.146E-01, 0.409E-05, 0.223E-01
     +, 0.827E+00, 0.534E-02, 0.145E+00, 0.880E-03
     +, 0.000E+00, 0.346E-01, 0.832E-01, 0.735E+00
     +, 0.809E-02, 0.139E-01, 0.124E+00, 0.133E-02, 0.0 /
      DATA (TMP(J, 67),J=1,ISHLS) /
     +  0.291E+00, 0.517E+00, 0.536E-01, 0.104E+00
     +, 0.118E-01, 0.229E-01, 0.224E-03, 0.139E-01
     +, 0.325E+00, 0.447E+00, 0.774E-01, 0.111E+00
     +, 0.104E-01, 0.144E-01, 0.465E-05, 0.222E-01
     +, 0.826E+00, 0.534E-02, 0.146E+00, 0.876E-03
     +, 0.000E+00, 0.349E-01, 0.831E-01, 0.734E+00
     +, 0.819E-02, 0.139E-01, 0.125E+00, 0.134E-02, 0.0 /
      DATA (TMP(J, 68),J=1,ISHLS) /
     +  0.291E+00, 0.516E+00, 0.538E-01, 0.104E+00
     +, 0.119E-01, 0.230E-01, 0.214E-03, 0.147E-01
     +, 0.327E+00, 0.445E+00, 0.780E-01, 0.111E+00
     +, 0.104E-01, 0.142E-01, 0.527E-05, 0.222E-01
     +, 0.826E+00, 0.534E-02, 0.146E+00, 0.871E-03
     +, 0.000E+00, 0.353E-01, 0.830E-01, 0.734E+00
     +, 0.830E-02, 0.139E-01, 0.125E+00, 0.134E-02, 0.0 /
      DATA (TMP(J, 69),J=1,ISHLS) /
     +  0.291E+00, 0.515E+00, 0.540E-01, 0.104E+00
     +, 0.120E-01, 0.232E-01, 0.205E-03, 0.155E-01
     +, 0.329E+00, 0.442E+00, 0.788E-01, 0.111E+00
     +, 0.104E-01, 0.140E-01, 0.598E-05, 0.222E-01
     +, 0.825E+00, 0.535E-02, 0.146E+00, 0.866E-03
     +, 0.000E+00, 0.357E-01, 0.830E-01, 0.733E+00
     +, 0.840E-02, 0.140E-01, 0.125E+00, 0.135E-02, 0.0 /
      DATA (TMP(J, 70),J=1,ISHLS) /
     +  0.292E+00, 0.514E+00, 0.542E-01, 0.105E+00
     +, 0.121E-01, 0.233E-01, 0.197E-03, 0.163E-01
     +, 0.331E+00, 0.439E+00, 0.794E-01, 0.110E+00
     +, 0.104E-01, 0.138E-01, 0.677E-05, 0.221E-01
     +, 0.825E+00, 0.536E-02, 0.147E+00, 0.861E-03
     +, 0.000E+00, 0.361E-01, 0.828E-01, 0.732E+00
     +, 0.851E-02, 0.140E-01, 0.125E+00, 0.136E-02, 0.0 /
      DATA (TMP(J, 71),J=1,ISHLS) /
     +  0.292E+00, 0.513E+00, 0.543E-01, 0.105E+00
     +, 0.121E-01, 0.235E-01, 0.189E-03, 0.172E-01
     +, 0.332E+00, 0.435E+00, 0.801E-01, 0.110E+00
     +, 0.110E-01, 0.145E-01, 0.765E-05, 0.221E-01
     +, 0.823E+00, 0.537E-02, 0.148E+00, 0.890E-03
     +, 0.649E-03, 0.365E-01, 0.827E-01, 0.730E+00
     +, 0.863E-02, 0.140E-01, 0.126E+00, 0.142E-02, 0.572E-03 /
      DATA (TMP(J, 72),J=1,ISHLS) /
     +  0.293E+00, 0.512E+00, 0.545E-01, 0.105E+00
     +, 0.122E-01, 0.237E-01, 0.181E-03, 0.181E-01
     +, 0.334E+00, 0.430E+00, 0.810E-01, 0.110E+00
     +, 0.116E-01, 0.152E-01, 0.858E-05, 0.220E-01
     +, 0.821E+00, 0.539E-02, 0.149E+00, 0.921E-03
     +, 0.160E-02, 0.369E-01, 0.825E-01, 0.728E+00
     +, 0.875E-02, 0.141E-01, 0.127E+00, 0.149E-02, 0.142E-02 /
      DATA (TMP(J, 73),J=1,ISHLS) /
     +  0.293E+00, 0.511E+00, 0.546E-01, 0.106E+00
     +, 0.124E-01, 0.239E-01, 0.174E-03, 0.191E-01
     +, 0.335E+00, 0.426E+00, 0.818E-01, 0.110E+00
     +, 0.122E-01, 0.158E-01, 0.968E-05, 0.220E-01
     +, 0.818E+00, 0.539E-02, 0.151E+00, 0.950E-03
     +, 0.278E-02, 0.373E-01, 0.822E-01, 0.726E+00
     +, 0.888E-02, 0.142E-01, 0.127E+00, 0.156E-02, 0.250E-02 /
      DATA (TMP(J, 74),J=1,ISHLS) /
     +  0.293E+00, 0.509E+00, 0.548E-01, 0.106E+00
     +, 0.125E-01, 0.242E-01, 0.167E-03, 0.202E-01
     +, 0.337E+00, 0.421E+00, 0.826E-01, 0.110E+00
     +, 0.128E-01, 0.165E-01, 0.109E-04, 0.219E-01
     +, 0.816E+00, 0.541E-02, 0.152E+00, 0.980E-03
     +, 0.420E-02, 0.377E-01, 0.820E-01, 0.723E+00
     +, 0.901E-02, 0.143E-01, 0.128E+00, 0.162E-02, 0.377E-02 /
      DATA (TMP(J, 75),J=1,ISHLS) /
     +  0.294E+00, 0.508E+00, 0.549E-01, 0.106E+00
     +, 0.126E-01, 0.244E-01, 0.161E-03, 0.213E-01
     +, 0.338E+00, 0.417E+00, 0.835E-01, 0.109E+00
     +, 0.134E-01, 0.170E-01, 0.122E-04, 0.219E-01
     +, 0.813E+00, 0.543E-02, 0.153E+00, 0.101E-02
     +, 0.581E-02, 0.381E-01, 0.817E-01, 0.721E+00
     +, 0.914E-02, 0.144E-01, 0.129E+00, 0.169E-02, 0.522E-02 /
      DATA (TMP(J, 76),J=1,ISHLS) /
     +  0.294E+00, 0.507E+00, 0.550E-01, 0.106E+00
     +, 0.127E-01, 0.246E-01, 0.155E-03, 0.225E-01
     +, 0.340E+00, 0.412E+00, 0.844E-01, 0.109E+00
     +, 0.140E-01, 0.176E-01, 0.137E-04, 0.218E-01
     +, 0.810E+00, 0.544E-02, 0.154E+00, 0.104E-02
     +, 0.762E-02, 0.384E-01, 0.814E-01, 0.718E+00
     +, 0.928E-02, 0.145E-01, 0.130E+00, 0.176E-02, 0.684E-02 /
      DATA (TMP(J, 77),J=1,ISHLS) /
     +  0.295E+00, 0.506E+00, 0.552E-01, 0.107E+00
     +, 0.128E-01, 0.249E-01, 0.149E-03, 0.239E-01
     +, 0.342E+00, 0.408E+00, 0.854E-01, 0.109E+00
     +, 0.144E-01, 0.178E-01, 0.153E-04, 0.217E-01
     +, 0.806E+00, 0.545E-02, 0.155E+00, 0.105E-02
     +, 0.109E-01, 0.388E-01, 0.810E-01, 0.714E+00
     +, 0.940E-02, 0.145E-01, 0.131E+00, 0.181E-02, 0.958E-02 /
      DATA (TMP(J, 78),J=1,ISHLS) /
     +  0.295E+00, 0.505E+00, 0.553E-01, 0.107E+00
     +, 0.129E-01, 0.251E-01, 0.144E-03, 0.253E-01
     +, 0.343E+00, 0.403E+00, 0.863E-01, 0.109E+00
     +, 0.151E-01, 0.184E-01, 0.172E-04, 0.217E-01
     +, 0.803E+00, 0.547E-02, 0.157E+00, 0.109E-02
     +, 0.125E-01, 0.393E-01, 0.807E-01, 0.711E+00
     +, 0.955E-02, 0.146E-01, 0.132E+00, 0.188E-02, 0.111E-01 /
      DATA (TMP(J, 79),J=1,ISHLS) /
     +  0.296E+00, 0.503E+00, 0.554E-01, 0.107E+00
     +, 0.130E-01, 0.254E-01, 0.140E-03, 0.270E-01
     +, 0.348E+00, 0.392E+00, 0.881E-01, 0.110E+00
     +, 0.158E-01, 0.190E-01, 0.192E-04, 0.216E-01
     +, 0.799E+00, 0.549E-02, 0.158E+00, 0.111E-02
     +, 0.148E-01, 0.397E-01, 0.803E-01, 0.708E+00
     +, 0.969E-02, 0.147E-01, 0.132E+00, 0.195E-02, 0.131E-01 /
      DATA (TMP(J, 80),J=1,ISHLS) /
     +  0.296E+00, 0.502E+00, 0.555E-01, 0.107E+00
     +, 0.131E-01, 0.256E-01, 0.133E-03, 0.283E-01
     +, 0.347E+00, 0.393E+00, 0.882E-01, 0.108E+00
     +, 0.163E-01, 0.194E-01, 0.214E-04, 0.216E-01
     +, 0.796E+00, 0.551E-02, 0.159E+00, 0.114E-02
     +, 0.165E-01, 0.401E-01, 0.801E-01, 0.705E+00
     +, 0.984E-02, 0.147E-01, 0.133E+00, 0.203E-02, 0.147E-01 /
      DATA (TMP(J, 81),J=1,ISHLS) /
     +  0.297E+00, 0.501E+00, 0.556E-01, 0.108E+00
     +, 0.132E-01, 0.259E-01, 0.129E-03, 0.301E-01
     +, 0.348E+00, 0.388E+00, 0.891E-01, 0.108E+00
     +, 0.170E-01, 0.200E-01, 0.240E-04, 0.215E-01
     +, 0.793E+00, 0.553E-02, 0.160E+00, 0.117E-02
     +, 0.183E-01, 0.406E-01, 0.798E-01, 0.703E+00
     +, 0.999E-02, 0.148E-01, 0.134E+00, 0.210E-02, 0.163E-01 /
      DATA (TMP(J, 82),J=1,ISHLS) /
     +  0.297E+00, 0.500E+00, 0.557E-01, 0.108E+00
     +, 0.133E-01, 0.261E-01, 0.124E-03, 0.319E-01
     +, 0.350E+00, 0.382E+00, 0.901E-01, 0.107E+00
     +, 0.176E-01, 0.206E-01, 0.267E-04, 0.215E-01
     +, 0.790E+00, 0.555E-02, 0.161E+00, 0.120E-02
     +, 0.200E-01, 0.411E-01, 0.795E-01, 0.700E+00
     +, 0.101E-01, 0.149E-01, 0.135E+00, 0.218E-02, 0.179E-01 /
      DATA (TMP(J, 83),J=1,ISHLS) /
     +  0.298E+00, 0.498E+00, 0.557E-01, 0.108E+00
     +, 0.134E-01, 0.263E-01, 0.120E-03, 0.338E-01
     +, 0.352E+00, 0.377E+00, 0.910E-01, 0.107E+00
     +, 0.183E-01, 0.212E-01, 0.298E-04, 0.215E-01
     +, 0.787E+00, 0.557E-02, 0.163E+00, 0.123E-02
     +, 0.216E-01, 0.415E-01, 0.792E-01, 0.697E+00
     +, 0.103E-01, 0.150E-01, 0.135E+00, 0.226E-02, 0.194E-01 /
      DATA (TMP(J, 84),J=1,ISHLS) /
     +  0.299E+00, 0.497E+00, 0.558E-01, 0.108E+00
     +, 0.135E-01, 0.266E-01, 0.116E-03, 0.358E-01
     +, 0.354E+00, 0.371E+00, 0.921E-01, 0.106E+00
     +, 0.190E-01, 0.218E-01, 0.332E-04, 0.214E-01
     +, 0.785E+00, 0.560E-02, 0.164E+00, 0.126E-02
     +, 0.232E-01, 0.420E-01, 0.789E-01, 0.694E+00
     +, 0.105E-01, 0.150E-01, 0.136E+00, 0.234E-02, 0.208E-01 /
      DATA (TMP(J, 85),J=1,ISHLS) /
     +  0.299E+00, 0.496E+00, 0.559E-01, 0.108E+00
     +, 0.135E-01, 0.268E-01, 0.112E-03, 0.381E-01
     +, 0.356E+00, 0.365E+00, 0.930E-01, 0.106E+00
     +, 0.197E-01, 0.223E-01, 0.369E-04, 0.214E-01
     +, 0.782E+00, 0.562E-02, 0.165E+00, 0.129E-02
     +, 0.247E-01, 0.426E-01, 0.786E-01, 0.692E+00
     +, 0.106E-01, 0.151E-01, 0.137E+00, 0.243E-02, 0.222E-01 /
      DATA (TMP(J, 86),J=1,ISHLS) /
     +  0.300E+00, 0.495E+00, 0.560E-01, 0.109E+00
     +, 0.136E-01, 0.270E-01, 0.108E-03, 0.405E-01
     +, 0.358E+00, 0.360E+00, 0.940E-01, 0.105E+00
     +, 0.204E-01, 0.227E-01, 0.411E-04, 0.214E-01
     +, 0.779E+00, 0.565E-02, 0.166E+00, 0.133E-02
     +, 0.262E-01, 0.431E-01, 0.783E-01, 0.689E+00
     +, 0.108E-01, 0.151E-01, 0.138E+00, 0.251E-02, 0.235E-01 /
      DATA (TMP(J, 87),J=1,ISHLS) /
     +  0.301E+00, 0.493E+00, 0.560E-01, 0.109E+00
     +, 0.137E-01, 0.273E-01, 0.105E-03, 0.430E-01
     +, 0.360E+00, 0.354E+00, 0.951E-01, 0.104E+00
     +, 0.211E-01, 0.232E-01, 0.457E-04, 0.214E-01
     +, 0.776E+00, 0.568E-02, 0.167E+00, 0.136E-02
     +, 0.276E-01, 0.437E-01, 0.781E-01, 0.687E+00
     +, 0.110E-01, 0.152E-01, 0.138E+00, 0.260E-02, 0.248E-01 /
      DATA (TMP(J, 88),J=1,ISHLS) /
     +  0.301E+00, 0.492E+00, 0.561E-01, 0.109E+00
     +, 0.138E-01, 0.275E-01, 0.101E-03, 0.456E-01
     +, 0.362E+00, 0.348E+00, 0.960E-01, 0.103E+00
     +, 0.218E-01, 0.236E-01, 0.508E-04, 0.214E-01
     +, 0.774E+00, 0.571E-02, 0.168E+00, 0.139E-02
     +, 0.291E-01, 0.425E-01, 0.779E-01, 0.685E+00
     +, 0.112E-01, 0.153E-01, 0.139E+00, 0.269E-02, 0.261E-01 /
      DATA (TMP(J, 89),J=1,ISHLS) /
     +  0.302E+00, 0.491E+00, 0.561E-01, 0.109E+00
     +, 0.139E-01, 0.277E-01, 0.978E-04, 0.485E-01
     +, 0.364E+00, 0.341E+00, 0.971E-01, 0.103E+00
     +, 0.224E-01, 0.240E-01, 0.564E-04, 0.215E-01
     +, 0.771E+00, 0.574E-02, 0.170E+00, 0.142E-02
     +, 0.305E-01, 0.448E-01, 0.775E-01, 0.681E+00
     +, 0.113E-01, 0.153E-01, 0.140E+00, 0.277E-02, 0.272E-01 /
      DATA (TMP(J, 90),J=1,ISHLS) /
     +  0.303E+00, 0.490E+00, 0.561E-01, 0.109E+00
     +, 0.139E-01, 0.280E-01, 0.946E-04, 0.515E-01
     +, 0.366E+00, 0.335E+00, 0.981E-01, 0.102E+00
     +, 0.231E-01, 0.243E-01, 0.626E-04, 0.215E-01
     +, 0.769E+00, 0.577E-02, 0.171E+00, 0.145E-02
     +, 0.319E-01, 0.454E-01, 0.773E-01, 0.679E+00
     +, 0.115E-01, 0.153E-01, 0.140E+00, 0.286E-02, 0.284E-01 /
      DATA (TMP(J, 91),J=1,ISHLS) /
     +  0.304E+00, 0.488E+00, 0.562E-01, 0.110E+00
     +, 0.140E-01, 0.282E-01, 0.913E-04, 0.547E-01
     +, 0.368E+00, 0.329E+00, 0.991E-01, 0.101E+00
     +, 0.237E-01, 0.244E-01, 0.695E-04, 0.215E-01
     +, 0.766E+00, 0.580E-02, 0.172E+00, 0.147E-02
     +, 0.326E-01, 0.460E-01, 0.770E-01, 0.677E+00
     +, 0.117E-01, 0.154E-01, 0.141E+00, 0.294E-02, 0.290E-01 /
      DATA (TMP(J, 92),J=1,ISHLS) /
     +  0.304E+00, 0.487E+00, 0.562E-01, 0.110E+00
     +, 0.140E-01, 0.284E-01, 0.882E-04, 0.582E-01
     +, 0.371E+00, 0.322E+00, 0.100E+00, 0.997E-01
     +, 0.243E-01, 0.245E-01, 0.772E-04, 0.215E-01
     +, 0.764E+00, 0.584E-02, 0.173E+00, 0.150E-02
     +, 0.335E-01, 0.466E-01, 0.768E-01, 0.675E+00
     +, 0.119E-01, 0.154E-01, 0.142E+00, 0.302E-02, 0.299E-01 /
      DATA (TMP(J, 93),J=1,ISHLS) /
     +  0.305E+00, 0.486E+00, 0.562E-01, 0.110E+00
     +, 0.141E-01, 0.286E-01, 0.852E-04, 0.619E-01
     +, 0.373E+00, 0.316E+00, 0.101E+00, 0.987E-01
     +, 0.249E-01, 0.246E-01, 0.855E-04, 0.216E-01
     +, 0.762E+00, 0.587E-02, 0.174E+00, 0.153E-02
     +, 0.344E-01, 0.473E-01, 0.766E-01, 0.673E+00
     +, 0.121E-01, 0.155E-01, 0.142E+00, 0.310E-02, 0.304E-01 /
      DATA (TMP(J, 94),J=1,ISHLS) /
     +  0.306E+00, 0.484E+00, 0.561E-01, 0.110E+00
     +, 0.141E-01, 0.298E-01, 0.822E-04, 0.658E-01
     +, 0.376E+00, 0.309E+00, 0.102E+00, 0.974E-01
     +, 0.254E-01, 0.246E-01, 0.949E-04, 0.216E-01
     +, 0.760E+00, 0.591E-02, 0.175E+00, 0.155E-02
     +, 0.351E-01, 0.480E-01, 0.764E-01, 0.671E+00
     +, 0.123E-01, 0.156E-01, 0.143E+00, 0.318E-02, 0.308E-01 /
      DATA (TMP(J, 95),J=1,ISHLS) /
     +  0.307E+00, 0.483E+00, 0.562E-01, 0.110E+00
     +, 0.142E-01, 0.291E-01, 0.793E-04, 0.700E-01
     +, 0.378E+00, 0.302E+00, 0.103E+00, 0.962E-01
     +, 0.259E-01, 0.245E-01, 0.105E-03, 0.216E-01
     +, 0.758E+00, 0.595E-02, 0.176E+00, 0.158E-02
     +, 0.359E-01, 0.486E-01, 0.763E-01, 0.669E+00
     +, 0.125E-01, 0.156E-01, 0.143E+00, 0.326E-02, 0.314E-01 /
      DATA (TMP(J, 96),J=1,ISHLS) /
     +  0.308E+00, 0.482E+00, 0.562E-01, 0.111E+00
     +, 0.142E-01, 0.293E-01, 0.764E-04, 0.745E-01
     +, 0.381E+00, 0.295E+00, 0.104E+00, 0.948E-01
     +, 0.265E-01, 0.245E-01, 0.116E-03, 0.217E-01
     +, 0.756E+00, 0.598E-02, 0.178E+00, 0.160E-02
     +, 0.369E-01, 0.493E-01, 0.761E-01, 0.667E+00
     +, 0.127E-01, 0.156E-01, 0.144E+00, 0.335E-02, 0.322E-01 /
      DATA (TMP(J, 97),J=1,ISHLS) /
     +  0.309E+00, 0.481E+00, 0.561E-01, 0.111E+00
     +, 0.143E-01, 0.295E-01, 0.736E-04, 0.792E-01
     +, 0.383E+00, 0.287E+00, 0.106E+00, 0.934E-01
     +, 0.271E-01, 0.244E-01, 0.129E-03, 0.218E-01
     +, 0.754E+00, 0.603E-02, 0.179E+00, 0.163E-02
     +, 0.376E-01, 0.500E-01, 0.757E-01, 0.665E+00
     +, 0.129E-01, 0.157E-01, 0.144E+00, 0.343E-02, 0.327E-01 /
      DATA (TMP(J, 98),J=1,ISHLS) /
     +  0.310E+00, 0.479E+00, 0.561E-01, 0.111E+00
     +, 0.143E-01, 0.297E-01, 0.708E-04, 0.843E-01
     +, 0.386E+00, 0.280E+00, 0.106E+00, 0.919E-01
     +, 0.276E-01, 0.242E-01, 0.143E-03, 0.218E-01
     +, 0.752E+00, 0.607E-02, 0.180E+00, 0.165E-02
     +, 0.381E-01, 0.507E-01, 0.756E-01, 0.664E+00
     +, 0.131E-01, 0.157E-01, 0.145E+00, 0.351E-02, 0.330E-01 /
      DATA (TMP(J, 99),J=1,ISHLS) /
     +  0.311E+00, 0.478E+00, 0.561E-01, 0.111E+00
     +, 0.143E-01, 0.299E-01, 0.679E-04, 0.897E-01
     +, 0.388E+00, 0.272E+00, 0.107E+00, 0.903E-01
     +, 0.281E-01, 0.240E-01, 0.158E-03, 0.219E-01
     +, 0.751E+00, 0.611E-02, 0.181E+00, 0.168E-02
     +, 0.388E-01, 0.515E-01, 0.754E-01, 0.662E+00
     +, 0.133E-01, 0.157E-01, 0.145E+00, 0.359E-02, 0.335E-01 /
      DATA (TMP(J,100),J=1,ISHLS) /
     +  0.312E+00, 0.477E+00, 0.560E-01, 0.111E+00
     +, 0.144E-01, 0.301E-01, 0.650E-04, 0.954E-01
     +, 0.391E+00, 0.264E+00, 0.108E+00, 0.885E-01
     +, 0.286E-01, 0.238E-01, 0.175E-03, 0.220E-01
     +, 0.749E+00, 0.616E-02, 0.182E+00, 0.170E-02
     +, 0.395E-01, 0.522E-01, 0.752E-01, 0.660E+00
     +, 0.135E-01, 0.158E-01, 0.146E+00, 0.367E-02, 0.339E-01 /
C
      DATA (TMP1(I,  1),I=1,4) / 4*0.0 /
      DATA (TMP1(I,  2),I=1,4) / 4*0.0 /
      DATA (TMP1(I,  3),I=1,4) / 4*0.0 /
      DATA (TMP1(I,  4),I=1,4) / 4*0.0 /
      DATA (TMP1(I,  5),I=1,4) /  0.170E-02, 3*0.0 /
      DATA (TMP1(I,  6),I=1,4) /  0.280E-02, 3*0.0 /
      DATA (TMP1(I,  7),I=1,4) /  0.520E-02, 3*0.0 /
      DATA (TMP1(I,  8),I=1,4) /  0.830E-02, 3*0.0 /
      DATA (TMP1(I,  9),I=1,4) /  0.130E-01, 3*0.0 /
      DATA (TMP1(I, 10),I=1,4) /  0.180E-01, 3*0.0 /
      DATA (TMP1(I, 11),I=1,4) /  0.230E-01, 3*0.0 /
      DATA (TMP1(I, 12),I=1,4) /
     +  0.300E-01, 0.290E-04, 0.120E-02, 0.120E-02 /
      DATA (TMP1(I, 13),I=1,4) /
     +  0.390E-01, 0.260E-04, 0.750E-03, 0.750E-03 /
      DATA (TMP1(I, 14),I=1,4) /
     +  0.500E-01, 0.300E-04, 0.370E-03, 0.380E-03 /
      DATA (TMP1(I, 15),I=1,4) /
     +  0.630E-01, 0.390E-04, 0.310E-03, 0.310E-03 /
      DATA (TMP1(I, 16),I=1,4) /
     +  0.780E-01, 0.740E-04, 0.260E-03, 0.260E-03 /
      DATA (TMP1(I, 17),I=1,4) /
     +  0.970E-01, 0.120E-03, 0.240E-03, 0.240E-03 /
      DATA (TMP1(I, 18),I=1,4) /
     +  0.118E+00, 0.180E-03, 0.220E-03, 0.220E-03 /
      DATA (TMP1(I, 19),I=1,4) /
     +  0.140E+00, 0.240E-03, 0.270E-03, 0.270E-03 /
      DATA (TMP1(I, 20),I=1,4) /
     +  0.163E+00, 0.310E-03, 0.330E-03, 0.330E-03 /
      DATA (TMP1(I, 21),I=1,4) /
     +  0.188E+00, 0.390E-03, 0.840E-03, 0.840E-03 /
      DATA (TMP1(I, 22),I=1,4) /
     +  0.214E+00, 0.470E-03, 0.150E-02, 0.150E-02 /
      DATA (TMP1(I, 23),I=1,4) /
     +  0.243E+00, 0.580E-03, 0.260E-02, 0.260E-02 /
      DATA (TMP1(I, 24),I=1,4) /
     +  0.275E+00, 0.710E-03, 0.370E-02, 0.370E-02 /
      DATA (TMP1(I, 25),I=1,4) /
     +  0.308E+00, 0.840E-03, 0.500E-02, 0.500E-02 /
      DATA (TMP1(I, 26),I=1,4) /
     +  0.340E+00, 0.100E-02, 0.630E-02, 0.630E-02 /
      DATA (TMP1(I, 27),I=1,4) /
     +  0.373E+00, 0.120E-02, 0.770E-02, 0.770E-02 /
      DATA (TMP1(I, 28),I=1,4) /
     +  0.406E+00, 0.140E-02, 0.860E-02, 0.930E-02 /
      DATA (TMP1(I, 29),I=1,4) /
     +  0.440E+00, 0.160E-02, 0.100E-01, 0.110E-01 /
      DATA (TMP1(I, 30),I=1,4) /
     +  0.474E+00, 0.180E-02, 0.110E-01, 0.120E-01 /
      DATA (TMP1(I, 31),I=1,4) /
     +  0.507E+00, 0.210E-02, 0.120E-01, 0.130E-01 /
      DATA (TMP1(I, 32),I=1,4) /
     +  0.535E+00, 0.240E-02, 0.130E-01, 0.150E-01 /
      DATA (TMP1(I, 33),I=1,4) /
     +  0.562E+00, 0.280E-02, 0.140E-01, 0.160E-01 /
      DATA (TMP1(I, 34),I=1,4) /
     +  0.589E+00, 0.320E-02, 0.160E-01, 0.180E-01 /
      DATA (TMP1(I, 35),I=1,4) /
     +  0.618E+00, 0.360E-02, 0.180E-01, 0.200E-01 /
      DATA (TMP1(I, 36),I=1,4) /
     +  0.643E+00, 0.410E-02, 0.200E-01, 0.220E-01 /
      DATA (TMP1(I, 37),I=1,4) /
     +  0.667E+00, 0.460E-02, 0.220E-01, 0.240E-01 /
      DATA (TMP1(I, 38),I=1,4) /
     +  0.690E+00, 0.510E-02, 0.240E-01, 0.260E-01 /
      DATA (TMP1(I, 39),I=1,4) /
     +  0.710E+00, 0.590E-02, 0.260E-01, 0.280E-01 /
      DATA (TMP1(I, 40),I=1,4) /
     +  0.730E+00, 0.680E-02, 0.280E-01, 0.310E-01 /
      DATA (TMP1(I, 41),I=1,4) /
     +  0.747E+00, 0.940E-02, 0.310E-01, 0.340E-01 /
      DATA (TMP1(I, 42),I=1,4) /
     +  0.765E+00, 0.100E-01, 0.340E-01, 0.370E-01 /
      DATA (TMP1(I, 43),I=1,4) /
     +  0.780E+00, 0.110E-01, 0.370E-01, 0.400E-01 /
      DATA (TMP1(I, 44),I=1,4) /
     +  0.794E+00, 0.120E-01, 0.400E-01, 0.430E-01 /
      DATA (TMP1(I, 45),I=1,4) /
     +  0.808E+00, 0.130E-01, 0.430E-01, 0.460E-01 /
      DATA (TMP1(I, 46),I=1,4) /
     +  0.820E+00, 0.140E-01, 0.470E-01, 0.490E-01 /
      DATA (TMP1(I, 47),I=1,4) /
     +  0.831E+00, 0.160E-01, 0.510E-01, 0.520E-01 /
      DATA (TMP1(I, 48),I=1,4) /
     +  0.843E+00, 0.180E-01, 0.560E-01, 0.560E-01 /
      DATA (TMP1(I, 49),I=1,4) /
     +  0.853E+00, 0.200E-01, 0.610E-01, 0.600E-01 /
      DATA (TMP1(I, 50),I=1,4) /
     +  0.862E+00, 0.370E-01, 0.650E-01, 0.640E-01 /
      DATA (TMP1(I, 51),I=1,4) /
     +  0.870E+00, 0.390E-01, 0.690E-01, 0.690E-01 /
      DATA (TMP1(I, 52),I=1,4) /
     +  0.877E+00, 0.410E-01, 0.740E-01, 0.740E-01 /
      DATA (TMP1(I, 53),I=1,4) /
     +  0.884E+00, 0.440E-01, 0.790E-01, 0.790E-01 /
      DATA (TMP1(I, 54),I=1,4) /
     +  0.891E+00, 0.460E-01, 0.830E-01, 0.850E-01 /
      DATA (TMP1(I, 55),I=1,4) /
     +  0.897E+00, 0.490E-01, 0.900E-01, 0.910E-01 /
      DATA (TMP1(I, 56),I=1,4) /
     +  0.902E+00, 0.520E-01, 0.960E-01, 0.970E-01 /
      DATA (TMP1(I, 57),I=1,4) /
     +  0.907E+00, 0.550E-01, 0.103E+00, 0.104E+00 /
      DATA (TMP1(I, 58),I=1,4) /
     +  0.912E+00, 0.580E-01, 0.110E+00, 0.111E+00 /
      DATA (TMP1(I, 59),I=1,4) /
     +  0.917E+00, 0.610E-01, 0.117E+00, 0.118E+00 /
      DATA (TMP1(I, 60),I=1,4) /
     +  0.921E+00, 0.640E-01, 0.124E+00, 0.125E+00 /
      DATA (TMP1(I, 61),I=1,4) /
     +  0.925E+00, 0.660E-01, 0.132E+00, 0.132E+00 /
      DATA (TMP1(I, 62),I=1,4) /
     +  0.929E+00, 0.710E-01, 0.140E+00, 0.139E+00 /
      DATA (TMP1(I, 63),I=1,4) /
     +  0.932E+00, 0.750E-01, 0.149E+00, 0.147E+00 /
      DATA (TMP1(I, 64),I=1,4) /
     +  0.935E+00, 0.790E-01, 0.158E+00, 0.155E+00 /
      DATA (TMP1(I, 65),I=1,4) /
     +  0.938E+00, 0.830E-01, 0.167E+00, 0.164E+00 /
      DATA (TMP1(I, 66),I=1,4) /
     +  0.941E+00, 0.890E-01, 0.178E+00, 0.174E+00 /
      DATA (TMP1(I, 67),I=1,4) /
     +  0.944E+00, 0.940E-01, 0.189E+00, 0.182E+00 /
      DATA (TMP1(I, 68),I=1,4) /
     +  0.947E+00, 0.100E+00, 0.200E+00, 0.192E+00 /
      DATA (TMP1(I, 69),I=1,4) /
     +  0.949E+00, 0.106E+00, 0.211E+00, 0.201E+00 /
      DATA (TMP1(I, 70),I=1,4) /
     +  0.951E+00, 0.112E+00, 0.222E+00, 0.210E+00 /
      DATA (TMP1(I, 71),I=1,4) /
     +  0.953E+00, 0.120E+00, 0.234E+00, 0.220E+00 /
      DATA (TMP1(I, 72),I=1,4) /
     +  0.955E+00, 0.128E+00, 0.246E+00, 0.231E+00 /
      DATA (TMP1(I, 73),I=1,4) /
     +  0.957E+00, 0.137E+00, 0.258E+00, 0.243E+00 /
      DATA (TMP1(I, 74),I=1,4) /
     +  0.958E+00, 0.147E+00, 0.270E+00, 0.255E+00 /
      DATA (TMP1(I, 75),I=1,4) /
     +  0.959E+00, 0.144E+00, 0.283E+00, 0.268E+00 /
      DATA (TMP1(I, 76),I=1,4) /
     +  0.961E+00, 0.130E+00, 0.295E+00, 0.281E+00 /
      DATA (TMP1(I, 77),I=1,4) /
     +  0.962E+00, 0.120E+00, 0.308E+00, 0.294E+00 /
      DATA (TMP1(I, 78),I=1,4) /
     +  0.963E+00, 0.114E+00, 0.321E+00, 0.306E+00 /
      DATA (TMP1(I, 79),I=1,4) /
     +  0.964E+00, 0.107E+00, 0.334E+00, 0.320E+00 /
      DATA (TMP1(I, 80),I=1,4) /
     +  0.965E+00, 0.107E+00, 0.347E+00, 0.333E+00 /
      DATA (TMP1(I, 81),I=1,4) /
     +  0.966E+00, 0.107E+00, 0.360E+00, 0.347E+00 /
      DATA (TMP1(I, 82),I=1,4) /
     +  0.967E+00, 0.112E+00, 0.373E+00, 0.360E+00 /
      DATA (TMP1(I, 83),I=1,4) /
     +  0.968E+00, 0.117E+00, 0.387E+00, 0.373E+00 /
      DATA (TMP1(I, 84),I=1,4) /
     +  0.968E+00, 0.122E+00, 0.401E+00, 0.386E+00 /
      DATA (TMP1(I, 85),I=1,4) /
     +  0.969E+00, 0.128E+00, 0.415E+00, 0.399E+00 /
      DATA (TMP1(I, 86),I=1,4) /
     +  0.969E+00, 0.134E+00, 0.429E+00, 0.411E+00 /
      DATA (TMP1(I, 87),I=1,4) /
     +  0.970E+00, 0.139E+00, 0.443E+00, 0.424E+00 /
      DATA (TMP1(I, 88),I=1,4) /
     +  0.970E+00, 0.146E+00, 0.456E+00, 0.437E+00 /
      DATA (TMP1(I, 89),I=1,4) /
     +  0.971E+00, 0.153E+00, 0.468E+00, 0.450E+00 /
      DATA (TMP1(I, 90),I=1,4) /
     +  0.971E+00, 0.161E+00, 0.479E+00, 0.463E+00 /
      DATA (TMP1(I, 91),I=1,4) /
     +  0.972E+00, 0.162E+00, 0.472E+00, 0.476E+00 /
      DATA (TMP1(I, 92),I=1,4) /
     +  0.972E+00, 0.176E+00, 0.467E+00, 0.489E+00 /
      DATA (TMP1(I, 93),I=1,4) /
     +  0.973E+00, 0.187E+00, 0.466E+00, 0.502E+00 /
      DATA (TMP1(I, 94),I=1,4) /
     +  0.973E+00, 0.205E+00, 0.464E+00, 0.514E+00 /
      DATA (TMP1(I, 95),I=1,4) /
     +  0.974E+00, 0.218E+00, 0.471E+00, 0.526E+00 /
      DATA (TMP1(I, 96),I=1,4) /
     +  0.974E+00, 0.228E+00, 0.479E+00, 0.539E+00 /
      DATA (TMP1(I, 97),I=1,4) /
     +  0.975E+00, 0.236E+00, 0.485E+00, 0.550E+00 /
      DATA (TMP1(I, 98),I=1,4) /
     +  0.975E+00, 0.244E+00, 0.490E+00, 0.560E+00 /
      DATA (TMP1(I, 99),I=1,4) /
     +  0.975E+00, 0.253E+00, 0.497E+00, 0.570E+00 /
      DATA (TMP1(I,100),I=1,4) /
     +  0.976E+00, 0.263E+00, 0.506E+00, 0.579E+00 /
C
      DO 20 J = 1,MAXELZ
         DO 10 I = 1,ISHLS
            GPHRAT(I,J) = TMP(I,J)
   10    CONTINUE
   20 CONTINUE
C
      DO 30 I = 1,29
         ISHLTR(I) = ITRN(I)
   30 CONTINUE
C
      DO 40 I = 1,MAXELZ
         GFLUPR(1,I) = TMP1(1,I)
         GFLUPR(2,I) = TMP1(2,I)
         GFLUPR(3,I) = TMP1(3,I)
         GFLUPR(4,I) = TMP1(4,I)
   40 CONTINUE
C
C K shell decay mode
      DO 50 I = 5,MAXELZ
         IF(I.LT.8) THEN
            IGNRFN(1,I) = 2
            IGNRFN(2,I) = 2
         ELSEIF(I.LT.81) THEN
            IGNRFN(1,I) = 3
            IGNRFN(2,I) = 4
         ELSE
            IGNRFN(1,I) = 3
            IGNRFN(2,I) = 3
         ENDIF
   50 CONTINUE
C L1 shell decay mode
      DO 60 I = 11,MAXELZ
         IGNRFN(3,I) = 4
         IF(I.LT.24) THEN
            IGNRFN(4,I) = 5
         ELSEIF(I.LT.49) THEN
            IGNRFN(4,I) = 9
         ELSEIF(I.LT.53) THEN
            IGNRFN(4,I) = 10
         ELSE
            IGNRFN(4,I) = 14
         ENDIF
   60 CONTINUE
C L2 shell decay mode
      DO 70 I = 12,MAXELZ
         IF(I.LT.13) THEN
            IGNRFN(5,I) = 5
            IGNRFN(6,I) = 5
         ELSEIF(I.LT.16) THEN
            IGNRFN(5,I) = 5
            IGNRFN(6,I) = 6
         ELSEIF(I.LT.29) THEN
            IGNRFN(5,I) = 6
            IGNRFN(6,I) = 7
         ELSE
            IGNRFN(5,I) = 8
            IGNRFN(6,I) = 9
         ENDIF
   70 CONTINUE
C L3 shell decay mode
      DO 80 I = 12,MAXELZ
         IF(I.LT.13) THEN
            IGNRFN(7,I) = 5
            IGNRFN(8,I) = 5
         ELSEIF(I.LT.16) THEN
            IGNRFN(7,I) = 5
            IGNRFN(8,I) = 7
         ELSEIF(I.LT.30) THEN
            IGNRFN(7,I) = 7
            IGNRFN(8,I) = 7
         ELSE
            IGNRFN(7,I) = 8
            IGNRFN(8,I) = 9
         ENDIF
   80 CONTINUE
      END
